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I Abstract. We present preliminary results for the curvature of the 

_^ pseudocritical line and susceptibilities in JV/ = 2 + 1 flavor QCD. The 

computations are carried out on lattice sizes of 16 3 x 4, at matching 

parameters of early work of the Bielefeld group. Emphasis is placed on 

d the control of systematic errors, by cross-validating results obtained by use 

of the Taylor expansion and measurements at imaginary chemical potential. 
O I To this end, we generalize the magnetic equation of state to the analysis of 

D the number density, and we extend it to imaginary values of the chemical 

i-Q potential. 
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1. Introduction 

Most likely, a large fraction of the participants in this meeting are familiar with 
the challenges faced by lattice QCD at non-zero baryon density [I]: the region 
of quark chemical potential /i q exceeding m^/2 remains prohibitively difficult, 
(T*j due to a severe sign problem [2]. 

For smaller chemical potential a variety of methods have been proposed to 
£> circumvent these difficulties, which have been quite successful in many respects. 

However, even in the easy region fi q < m n /2 there are challenges to face: for 
instance, it has been noted that most current estimates of the QCD endpoint 
lie dangerously close to the 'forbidden region' \x q > m 7T /2; the qualitative 
agreement among different results on the slope of the critical line still needs 
to be corroborated by some convincing continuum estimates obtained with 
different methods. And, there is a growing interest in the values of higher 
order susceptibilities, particularly at larger values of the chemical potential, 
while the current results are limited to the lowest order ones 13). So there is 
room for improvement, which can perhaps be achieved by mixing and matching 
different approaches. 
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Figure 1. A sketchy view of the phase diagram of QCD for two massless flavors 
in the temperature, imaginary chemical potential plane, with the superimposed 
trajectories of our simulations. The blue points indicate the Roberge- Weiss 
endpoint (left). The chiral condensate measured along such trajectories, 
plotted against the rescaled chemical potential [i 2 — ^(T) (right). Note the 
(approximate) universal behaviour of the condensate. 



2. Two complemetary approaches: Taylor expansion and imaginary 
chemical potential 

One way to handle the physics of small chemical potential relies on Taylor 
expansion |4|. This approach is well understood mathematically, however it 
only works inside its convergency circle, whose radius is unknown. Moreover 
even inside the circle it is hard to judge when convergence has been achieved. 
Alternatively, analytic continuation from imaginary chemical potential - where 
conventional simulations are possible - to real chemical potential can in 
principle be extended to the entire analytic domain, but the danger here is 
a possible lack of control over the continuation itself: for instance, terms which 
are subleading for imaginary chemical potential might well become leading 
for real chemical potential. The strengths of the two methods are somewhat 
complementary, and their weaknesses different: this suggests to develop a 
hybrid method which combines the two, and this note reports on our ongoing 
effort in this direction J5j. 

3. Results 

Our study for staggered 2+1 flavors on lattices of size 16 3 x 4 clearly is 
exploratory. However, the quark mass values are close to the ones realized 
in nature. The Goldstone pion mass is tuned to about 220 MeV and the kaon 
acquires its physical mass. This corresponds to a (degenerate) light to strange 
quark mass ratio of mi/m a = 1/10. The action utilized is the p4fat3 action 
for which Taylor coefficients computed at \i = are available in our range of 
temperatures and at our quark masses |4|. The computations were carried out 
at three temperatures above the critical temperature for two massless flavors for 
the considered lattice spacings: T = 205, 210, 218 MeV. In Fig.l, left we sketch 
the simulations' trajectories, alongside the critical line in the chiral limit, in the 
Hi,T plane. The sketch evidences the /// periodicity as well as the symmetry 
under reflection fj,j — > —fJ>i; the actual simulations were performed only at 
positive three- flavor degenerate ///, and within the first half-period. The blue 
points mark the endpoint of the Roberge Weiss phase transition. The light 
quark chiral condensate, plotted in Fig.l, right, is showing a fairly rapid rise 
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Figure 2. The fits motivated by the magnetic Equation of State for the chiral 
condensate,superimposed to the data point. 



when fij is increased, signaling the pseudocritical point for the transition to the 
broken phase (marked SU(2)y in the diagram). We plot the results against a 
rescaled chemical potential to confirm the feasibility of a scaling study based 
on the magnetic Equation of State. 

3.1. The magnetic Equation of State 
Consider now the magnetic equation of state 



M = h 1/s f G (z) z = thW 



(1) 



where M and h are identified with the light chiral condensate and the light 
bare mass: M = m s ("tpijj)i, h = j^-^- respectively, fo is a scaling function 
which only depends on the universality class of the transition, and t = defines 
the critical line. Our light mass does not change: if scaling (i.e. Eq. 1) holds 
true the chiral condensate results plotted versus jj, 2 — [J? C (T), where ^ 2 {T) is the 
critical chemical potential for each temperature T, should fall on a single curve. 
In Fig.l we see that this is indeed happening, demonstrating the universal 
behaviour of the chiral condensate at least in the range of parameters we are 
exploring. fi 2 (T) is found by searching for the optimal match between results 
for different temperatures, and compares favorably with alternative estimates. 



3.2. Curvature of the critical line from the magnetic EoS 
For small chemical potential t can be explicitely written as 
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where T c is the critical temperature in the chiral limit. The standard condition 
/g(0) = /g(0) = 1 fi xes the normalization to an d ho which have been 
determined in Rcf. |6|. As a working hypothesis we assume that to and /io 
do not depend on the temperature, and we check a posteriori that this indeed 
is the case, k - the curvature of the critical line - can then be computed 
from the Equation of State, using the explicit forms available for the two 
plausible universality classes 0{2) and 0(4), and the results contrasted with 
those obtained in Ref. [7]. 

From the EoS for the order parameter one can derive analogous relations 
for the magnetic susceptibility, as well as for the number density n and its 
susceptibility |8|, leading to alternative estimates of k. To fix the notation, let 



Table 1. The results for the curvature of the critical line k (Eq.(2)) from the 
fits to the magnetic EoS for the chiral condensate (k^), and from the fits to 
the number density (k n ). 
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Figure 3. The fits motivated by the magnetic Equation of State for c\ 
and C2 (related to quark number and its susceptibility respectively, see text) 
superimposed to the data points. 

us write the Taylor expansion of the pressure in terms of the quark chemical 
potentials around a generic jXq = (fi u ,o, Hdfi-, Ms,o) for our Nj = 2 + 1 system: 
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with the abbreviation jj,„ = \i q jT and the coefficients 



Qk Ql gn 



Cklr 



kW.nl dfl u k dfi d l dfi 



r n 
s 



p_ 



(3) 



evaluated at uq. We will in particular discuss c\ = cioo and ci = C200, which 
are related to the particle number and its susceptibility. In Fig. 2 and Fig. 3 we 
show the results of the EoS fits and in Table 1 we record the k values returned 
by the fits for the chiral condensate and the quark number (the coefficient c\). 
While the fit results, at least for 0(2), do not seem to be inconsistent, it is clear 
that at values for fij close to the Roberge- Weiss limit other scaling studies need 
to be performed, work for the future. 



3.3. The Taylor expansion 

How far in \i we can trust the Taylor expansion? A direct comparison of the 

number density of e.g. the u quark given by the partial sums of its Taylor series: 

"KA) = 5Z kc kin{i*u ~ /vo) fc_1 (fJ-d ~ iJ-dfi) l (ft-s ~ As,o) n - ( 4 ) 

k,l,n 







Figure 4. The results for the number density at imaginary chemical potential, 
contrasted with the partial sums of the Taylor series of the indicated order. 



with results obtained directly at imaginary fi might give the answer. We show 
the results for increasing temperatures in Fig. 4: for T = 205 MeV a good 
convergence is achieved indicating a small contribution from the 8 order terms 
and superior ones. At the highest temperature instead the results suggest 
that a 6 order expansion is not capable to describe the data at [ij ^0.2. 
The simultaneous consideration of Taylor expansion and imaginary \x results 
naturally suggests a strategy to extract higher order terms of the expansions - 
or, equivalently, higher order susceptibilities - withouth having to implement 
the cumbersome higher order derivatives. At the same time one can place lower 
bounds on the radius of convergence simply by assessing by eye the extent of 
the region where actual convergence is achieved. 

4. Parting comments 

We have presented an overview of results obtained by performing simulations 
at zero and non-zero imaginary chemical potential. We have noted that an 
analysis based on the magnetic equation of state for the chiral condensate, the 
number density and their susceptibilities at imaginary chemical potential gives 
information on the universality class of the transition for two massless flavors, 
and on the slope of the chiral line. We have argued that the simultaneous 
consideration of results from the Taylor expansion at /i = and at imaginary 
/i can place a lower bound on the radius of convergence at /x = 0, and facilitate 
the computation of higher order coefficients. 
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